ZU-TH 9/98 
gr-qc 9805065 



OO 
On 



O 

in 



X 



THE GENERALIZED THIN-SANDWICH PROBLEM 
AND ITS LOCAL SOLVABILITY 

DOMENICO GlULINI 



Institut fur theoretische Physik, Universitat Zurich 
Winterthurerstrasse 190, CH-8057 Zurich, Switzerland 
e-mail: giulini@physik.unizh.ch 

00 



^ ' Abstract 

in 



We consider Einstein Gravity coupled to matter consisting of a gauge field with 
any compact gauge group and coupled scalar fields. We investigate under what 
■ conditions a free specification of a spatial field configuration and its time derivative 

qq ' determines a solution to the field equations (thin-sandwich problem). We establish 

, sufficient conditions under which the thin-sandwich problem can be solved locally 

""^ ■ in field space. 

cr 



Introduction 

In this paper we consider the initial value problem for Einstein gravity plus matter 
in spacetimes ExR, where £ is a closed orientable 3-manifold. We are interested 
in the question of how to find initial data which satisfy the constraints. The most 
popular approach here is a powerful method devised by Lichnerowicz, Choquet- 
Bruhat, York and others, henceforth referred to as the "conformal method". (See 
[I] for a brief review and [CY] for more details.) Of the gravitational variables it 
allows to freely specify the conformal class of the initial 3-metric, the conformally 
rescaled transverse-traceless components of the extrinsic curvature and a constant 
trace thereof (i.e. £ must have constant mean curvature). Given these data, the 
constraints turn into a quasilinear elliptic system of second order for the confor- 
mal factor (scalar function) and the transverse momentum (vector field), which 
decouples due to the constant mean- curvature condition. The disadvantages of this 
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method are that it does not easily generalize to data of variable mean curvature and 
that it does not allow to control the local scales of the physical quantities initially, 
since the freely specifiable data (gravitational and non-gravitational) are related to 
the actual physical quantities by some rescalings with suitable powers of the confor- 
mal factor. In particular, one has no control over the conformal part of the initial 
3-geometry. 

In this paper we are concerned with the so called "thin-sandwich method", 
which differs from the one just mentioned insofar as it aims to define solutions 
to the Einstein equations by a free specification of the initial field configuration 
and its coordinate-time-derivative. The constraints are now read as equations for 
the gauge parameters (lapse, shift, ..). From the conformal point of view this 
means that one tries to trade in the freedom to specify the gauge parameters for 
the freedom to specify the conformal part of the metric and the longitudinal part 
of the momentum. The disadvantages mentioned above would then be overcome, 
but unfortunately the equations (for the gauge parameters) turn out to be non- 
elliptic^ in general [BO][CF]. However, for certain open subsets of initial data 
they are elliptic and can be locally solved. This was first shown in [BF] and will be 
shown in a wider context with dynamical matter here. 

We note that historically this approach arose from the question (then formu- 
lated as a conjecture) of whether the specification of two 3-geometries uniquely de- 
termine an interpolating Einstein space-time (thick-sandwich problem) . For nearby 
geometries infinitesimally close in time this turns into the thin-sandwich problem 
(see [W] (chapter 4) and [BSW]) which we now describe in more detail. 

In a space-time neighborhood E x R of the Cauchy surface E, we use the 
standard parametrization of the space-time metric 

g M = - a 2 dt®dt + g ab (dx a + f3 a dt) <g> (dx b + f3 b dt), (1.1) 

where g is the (t-dependent) Riemannian metric of E and a, (3 are (t-dependent) 
scalar- and vector fields on E, known as lapse and shift. The extrinsic curvature 
reads 

K=±{d t -L f j)g, (1.2) 

By ellipticity of non-linear differential operators one means the ellipticity of its linearization, 
which depends on the point (in field space) about which one linearizes. The usual statement 
that the thin-sandwich equations are not elliptic merely asserts the existence of points where the 
linearization is not elliptic, but not that the domain of ellipticity is empty. 
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where Lp denotes the Lie-derivative along (3. Written in terms of (g, K), the con- 
straints do not depend on a and (3 and hence constrain the set of allowed values for 
the data (g, K). The constraints read 



where T is the energy-momentum tensor of the matter and _L denotes the component 
along the future pointing normal n of E. 

Alternatively, one may write the constraints in terms of g and g := dtg by 
replacing K via (1.2). Then they explicitly involve a and (3 and one may ask 
whether it is possible to freely specify (g, g) and let the constraints determine a 
and f3. This is the thin-sandwich problem (TSP). If we abbreviate ^> := (<?,<?), 
X := (a, f3), the constraints take the form of the thin-sandwich equation (TSE): 



The TSP now asks for existence and uniqueness of solutions of the TSE, read as 
equation for X given Once this is solved, we can construct (g, K) satisfying (1.3- 
4), which uniquely determine space-time via the Einstein evolution equations [CY]. 

It has long been discussed in the literature that, in general, existence and 
uniqueness should fail [BO], although some arguments merely showed this under 
the additional assumption (a priori) of constant lapse function a [CF] (see also [G] 
for a related issue). However, more recently it was shown that given a solution 
(*, X) of the TSE which satisfies certain bounds on geometric quantities and which 
admits no nontrivial solutions of the spatially projected Killing equation, there exist 
unique solutions X'(^l') for all in a neighbourhood of ^> . This was achieved by 
an implicit-function-theorem for a reduced version of (1.5) with already eliminated 
lapse function. 

Note that in this formulation of the TSP the right hand sides of the constraints, 
that is Tj__i_ and , are assumed given. These are not the components of T that an 
observer along d t would measure. The relation between the two sets of components 
involve a and (3. For example, if the matter is represented by some dynamical 
field (f>, the quantities T±± and T" cannot be calculated from the initial data ((f), (f>) 
without the use of a and (3. Hence there is a certain inconsistency in the traditional 



K ab K ab -(K*) 2 -R=-2T ±± , 
V\KZ-8 a b Kt)=Ti, 



(1.3) 
(1.4) 



F[V,X] = 0. 



(1.5) 
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formulation of the TSP, in that it eliminates any appearance of the normal n in 
favour of dt, lapse and shift on the left side (the gravitational part), but not on the 
right side (matter part) of the constraints. This inconsistency was already felt by 
others (see e.g. section IV of [BO]), but no alternative formulation was hitherto 
attempted. 

In this paper we consider a generalized thin-sandwich problem (GTSP) for 
the initial data of the full system of coupled gravitational and matter fields, which 
avoids the difficulty just mentioned. As matter we shall consider coupled systems 
of scalar and gauge fields, further specified below. By $ we shall collectively denote 
all dynamical fields of the theory. We ask: Under what conditions do input data 
\1/ := ($, $) uniquely specify a solution to the Einstein-matter equations? In the 
same fashion as for (1.5), one obtains a now generalized thin-sandwich equation 
(GTSE) from which one tries to determine the "gauge parameters" X given the 
data One non-trivial aspect of our generalization is due to the possible presence 
of gauge matter fields. In this case X comprises lapse, shift and additional functions 
with values in the Lie-algebra of the gauge group. Our main result will consist of 
an implicit-function-theorem for this extended set of variables, which for our GTSP 
is precisely analogous to the result proven in [BF] for the traditional form of the 
thin-sandwich problem. But note that the two formulations differ even without 
gauge fields. 



The Generalized Framework 

The dynamical fields we consider involve the gravitational field, a gauge field with 
compact gauge group G of dimension N and an M-component scalar field with 
values in an associated R M -vector-bundle. It couples to both previous fields in the 
standard minimal fashion. For simplicity we assume the G-principal bundle to be 
trivial. Since the frame bundle of any orientable 3-manifold is always trivial, we 
may once and for all choose global trivialisations of these bundles and represent 
fields by their globally defined component-fields on space-time. For fixed time, a 
configuration of fields is given by the 6 + 3iV + M component-fields on E, 

$ A :=G7a6,^,<n, (2.1) 

where indices fx, v.. denote components in the Lie-algebra LG and a, (3.. denote 



4 



components in R . Hence we think of a field configuration as mapping 

$ : E — ► GL(3, R)/SO(3) x K 3N x R M , (2.2) 

where we identify the space of symmetric, positive definite matrices, in which g a b 
is valued, with the first factor space on the right hand side. The total target space, 
whose dimension is 6 + 3iV + M, will be denoted by 0, and the space of mappings 
E — > (to be further specified) by M. 

Compactness of the gauge-group G implies the existence of G-invariant, sym- 
metric, positive definite, bilinear forms and h a p on LG and R M respectively. 
The class of models we shall consider here are characterized by the Lagrange four- 
form 

C=\* - | A * n " ~ lh af, V<l> a A *V(f)p - W, (2.3) 

where * is the Hodge-duality map wrt. and O is the curvature of A. For nota- 
tional simplicity we shall denote all the covariant derivatives acting on sections in 
the various vector bundles by the same symbol V. In general it will therefore involve 
the Christoffel symbols of g as well as the gauge connection A in the appropriate 
representation of LG. The potential W depends on the fields and its first spatial 
derivatives. Its precise form is not important, except that we need to explicitly 
exclude 2— (or higher) derivative couplings, in particular, the so called conformal 
coupling of the scalar and gravitational field. The reason for this will be explained 
below. The Hamiltonian constraint for (2.3) has the general form 

H = \G AB V A V B + U = 0, (2.4) 

where the potential is the following sum: 

U = -R+ \k^g ac g bd W ai p? cd + h a pg ab V ' a a V \<\P + W, (2.5) 

whose terms represent the contributions of the gravitational, gauge, and scalar 
fields (two terms) respectively to the potential energy. R denotes the Ricci-scalar 
for g. The "canonical velocities", V A , can be written in terms of the "coordinate 
velocities" <b A := d t & A and the "gauge-parameters" a and £ := (/?, A). The general 
structure is (compare (1.2)) 

v = ±r = ±(<i> + fz), (2.6) 
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where /e represents the motion generated by the infinitesimal diffeomorphism and 
gauge-transformation with parameters (3 and A respectively. Resolved in terms of 
the individual fields (2.1), the components of T, and hence of fe, read: 

Tab = {jab ~ 2V( a /3 6 ) , (2.7) 

K = K- P b nl - V a A^ , (2.8) 
r« = 0« - p a v a <p a + A"p°,^, (2.9) 

where p denotes the representation of LG in gl(M, R). 

Finally, following the decomposition of 6 as cartesian product, the "kinetic- 
energy-metric" Gab on 6 which appears in (2.4) has the following block-structure: 

G AB = G abcd © k^g ab © h af) , (2.10) 

where the first 6x6 block is given by the DeWitt-metric 

G abcd = 1 ( g ac g bd + g ad g bc _ 2 g ab g cd ) , (2.11) 

which is a Lorentz metric of signature (1,5). Hence Gab itself is a Lorentz metric of 
signature (1, 5 + 3N + M) on the manifold 0, which is homeomorphic to R_ 6 + 3Ar + M . 
We shall sometimes denote this metric simply by G and write G(-, •) for the inner 
product. We will see that the Lorentzian signature of G is the important feature 
on which the proofs of our main results rest. This is also the reason why we had to 
exclude higher derivative (e.g. conformal) couplings of the scalar and gravitational 
fields, since they will in general destroy this signature structure [K]. On the other 
hand, our proofs will still apply to more complicated self-couplings of the scalar field. 
For example, non-linear a-models would be allowed, since here the target space 
metric of the scalar field, h a p, simply becomes ^"-dependent, which is unimportant 
to our proofs as long as it stays positive definite. 

The (undensitized) momenta of the field $ A are just given by the covariant 
components - with respect to G - of the velocities: Pa := GabV b . For the indi- 
vidual fields we write Pa = {n ah ', 7r^, 7r a ). In the canonical theory, the phase space 
function that generates infinitesimal diffeomorphisms and gauge-transformations 
with parameter £' is given by 

V? := [ dfiP A f^, (2.12) 
Jt. 
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where here and below we set dji := ^det{g a b}d 3 x for the measure on E. For 
completeness we remark that, using (2.7-9), a straightforward calculation yields the 
following familiar Poisson-bracket relation, which involves the curvature tensor Q 
of the gauge field on E: 

{V e ,Vt»} = P e », (2.13) 

where f" = (/?'", A'") reads 

p"' = -[p',p"}, (2.14) 
\'" = [\>,\"}-n(f3',f3"). (2.15) 



The diffeomorphism- and Gaufi-constraints are just given by V^> = 0V£', which 
may be expressed by saying that the velocity field V is L 2 G-orthogonal to all "ver- 
tical" vector fields f^>. Hence we must have 

0= / P A f£dn= j ±G(T,f e )d»=: j (g ab (3' a V b + k^X'^dfi, (2.16) 

for all, say C°°, vector fields and LG-valued functions A'. This is equivalent to 
T> a = (diffeomorphism constraint) and = (Gaufi constraint). Explicitly we 
get 

v* = 2G abcd v b (^r cd ) - ^ a V^ cd r^ - ^ a6 (v 6 a )/i Q/3 r^ , (2.17) 
= v a (£</°*r£) + ihrfiptPh^ri . (2.18) 

Given $ and $ we now have the 4 + iV equations = = V a = for the 
4 + AT unknowns ct,/3 a , A M . The first step consists of inserting (2.6) into (2.4) and 
solving for a 2 : 



For this to make sense the right-hand side must be positive. But for the following 
analysis it turns out that we need to put the following stronger 



Condition 1 (a priori): 



U>0, (2.20) 
G(r,T)<0. (2.21) 
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Note that (2.20) just involves the initial data, whereas (2.21) contains as well the 
unknowns (3 and A (hence "a priori"). Note that (2.21) says that the system must, 
at each point of S, move in a "timelike" direction with respect to the Lorentz 
metric G. The need for such an a priori bound implies that our results will only 
be perturbative. Given these bounds, we set a equal to the positive square root 
of the right hand side of (2.19). We can then eliminate a from (2. 17), (2. 18) and 
obtain a set of 3 + N equations for the 3 + N unknowns £ = (/?, A), which we call 
the generalized reduced thin-sandwich equation (GRTSE): 

F[*,£] = 0. (2.22) 



Now consider the following functional over configurations satisfying Condi- 
tion 1: 

S[*,£]= [ y/-2UG(T,T) dp, (2.23) 

JT, 

then solutions to the GRTSE are stationary points with respect to variations in £. 
More precisely, let D 2 denote the partial (functional-) derivative in the second argu- 
ment and denote the L 2 inner product of £i and £2 by (C1IC2) : = / s dfi(k^ lv ,Xi fl X2 1 ' + 
9abPt$)- We then have 



Lemma 1. 



D 2 S[*,Z](0 = -(e\F[*,Z}} 



(2.24) 



Proof. For s e (-e, e) set t/(s) = £ + s£' and T(s) = <& + f v ( s ). Hence ^| s=0 r(s) 
fa. Then, recalling (2.16), 



d_ 

ds 



c(r,/ e )^ = -(e / |F[*,e]) 



(2.25) 
□ 
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Main Results 



In this section we are mainly concerned with the linearization of the GRTSE, ex- 
cept at the end where we will discuss global uniqueness. The corresponding linear 
operator will be called L, without explicit indication that it depends on \1/ and £. 
It is defined by 



as s=o Jy, 



d 
ds 



s=0 



G(fe,r( 8 )) 

a(s) 



(3.1) 



where T(s) = $ + / c+s? » and a(s) = [-G(T(s),T(s))/2U] 1 / 2 . Setting T(s = 0) = T, 
a(s = 0) = a, and noting that ^t\ s=0 ^( s ) = /$"> we S e t 



d 
ds 



a(s) 



G(ft»,T) G(/ r ,T) 



<:i 2atf ' G(r,r) 



where we used (2.4) to eliminate U in the last step. Hence 

G(/ € ',r)G(/ € »,r) 



ds 
where 



s=0 



g(A',r( 8 )) i 

a(s) a 



G{fc,h»)- 



G(r,r) 
G(/ e ,r) 



(3.2) 



(3.3) 



(3.4) 



G(r,r) 

is the G-orthogonal projection of perpendicular to V. This leads to the following 
expression for L's matrix elements: 



W) = J^iG(f^f^), 



(3.5) 



where a is the square-root of the rhs. of (2.19). If (2.21) holds, V is "timelike" 
(pointwise on E) and hence the / s are "spacelike" or zero. Since the metric G is 
Lorentzian, it is positive definite on "spacelike" vectors. Hence we have shown 

Lemma 2. Suppose Condition 1 holds, then L is self-adjoint and non-negative. 
Furthermore, £' G kernel (L) 3 k : £ — > R such that 



(3.6) 



Remark: Symmetry is expected, for consider the rhs. of (2.23) as functional of \& 
and r, denoting it by S[$, T], then the calculation of the rhs. of (3.5) was just that 

of -d$s[v, r](r', r") for r = / e and r" = fo. 

We will next show that under the same hypotheses L is in fact elliptic. For 
this we will need 
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Lemma 3. If Condition 1 holds, 7t~ ab is a positive or negative definite matrix. 

PROOF. Condition 1 implies G AB P A P B < =>• G abcd n ab n cd = 2(7r a& 7r a6 -±«) 2 ) < 
0. Choosing a frame where g ab = 5 ab and n ab = diag (pi,P2,P3), this is equivalent 
to the following condition on the eigenvalue-vector p := (pi,P2,P3)' 



n ■ p 



> 



(3.7) 



where n = (1, 1, l)/y/3, which means that p lies in the interior of the double cone 
with axis along n and opening angle 9 < cos _1 (2/3) 1/ ' 2 about its axis. This cone 
just touches the walls of the positive and negative octants along the bisecting lines, 
and since p must be in its interior, all eigenvalues are either strictly positive or 
strictly negative. □ 

Proposition 1. Suppose Condition 1 holds, then the second order differential op- 
erator L is elliptic. 

Proof. We shall calculate the principal symbol of L. For this we need to go back to 
the explicit formulae (2.17) and (2.18) for the full non-linear problem and explicitly 
linearize them, but keeping track only of the highest (second) derivatives. By = we 
shall denote equality in the second derivative terms. We set again T(s) = $ + /^ +s ^ 
etc., with £' = (/?', A') . It will be convenient to express things in terms of the 
momenta, using olPa = Gab^ B , and accordingly write (3.2) in the form 



Then 



d 
ds 



s=0 



V a = d b 



d 
ds 



s=0 



a(s) = 



PaH 
2U 



(3.8) 



aU 



-7T 



G abcd d c f3' d 



a 



1 

a 



2TT ac n de g db d c d e l3'° + n ab n c „d b d c \ 
U 



+ g bc d b d c (3 ,a - g ac d c d b (3' c 



(3.9) 



d 
ds 



s=0 



1 

a 



^lk^K--9 ah db\" 

2aU a 

ib 



2k^^g bc 7i cd d a d d /3' b + k^y b u d a d b X' v 
2U 



+ g ab d a d b X 



(3.10) 
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Replacing d a — > k a we can just read-off the matrix of the principal symbol cr(k) in 
the general form 

r 

(3.11) 



< 


< 


< 


< 



where we have chosen to order the N + 3 rows and columns so that we first count 
the iV components of A' and then the 3 components of (3' . In order to calculate the 
determinant we make the following simplifications: We call Ti ah kb 



--: p a , 7v^k a =: tt m , 



7r M := k^ v -K v and choose a spatial frame where g a b = 8 a b, k a = 
P a = (Pi)P2,0). Then (3.11) reads explicitly 



|fc||,0,0) and 



a(k) = -- 







7T M pi 

u 


7T M p 2 

u 





1 


u 


2p\ 
U 


2piP2 





a 


U 


2piP2 

U 


2 4 + \\k\\ 2 
















\\k\\ 2 



(3.12) 



Now, for k 7^ 0, pi = \\k\\ir(k, k) : where k := k/\\k\\. Lemma 3 then ensures that 
Pi 7^ 0. In order to calculate det{cr(/c)}, we simplify this matrix as follows: We 
subtract times the N + 1— row from the fx— row, for each 1 < \i < N, and 

also subtract — times the N + 1— row from the N + 2— row. The resulting matrix 

pi ° 

reads 





~\\k\\ 2 5» 











1 




2p? 


2piP2 







u 


U 




a 








\\k\\ 2 
















\\k\\ 2 



(3.13) 



Its determinant, which equals that of cr(k), is now easily calculated: 



det{a(/c)} 



2pf 
U 



2(N+2) 



r 


N+3 




= 2 


[ ll^ll 2 ] 


a 


a 



N+3 



U 



(3.14) 



Lemma 3 implies that this is zero k = 0, which finally proves ellipticity of L. □ 

The results obtained so far suffice to deduce an implicit-function-theorem. To 
state it precisely, we need to choose appropriate function spaces. It is natural 
to choose Sobolev spaces since they are also used to show existence for the time 
evolution [CY]. To begin with, it is convenient to summarize the order of spatial 
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differentiation by which the various fields enter the quantities V, U, a and hence 
the GRTSE, by the following matrix: 





r 

J. 


77 

\J 


LX 


ORTSE 


Hab 


i 

J. 


9 


9 


Q 
o 


a 


1 


1 


1 


2 


(f) a 


1 


1 


1 


2 


dab 










1 












1 












1 


/3 a 


1 




1 


2 


\^ 


1 




1 


2 



(3.15) 



Note that we assumed that U only contained first derivatives of the matter-fields, 
whereas it contains second derivatives of the gravitational field through the Ricci- 
scalar. Hence g a b enters the GRTSE thrice differentiated. ^ 

By H n (V) we shall denote the Sobolev space of V- valued functions on E with 
L 2 -norm in the first n derivatives (i.e. generalizing H n = W n ' 2 using an inner 
product in V). We shall have V = T 2 ° for g ab and g ah , V = Tf <g> LG for A^ and A%, 
V = R M for <p a and <p a , V = Tq for f3 a and V = LG for A M . The inner products for 
the various V's are just as in the metric G, except for the gravitational field where 
instead of G abcd , which is not positive definite, we choose the positive definite form 
gacgbd ( com p are (2.11)). Now we define the Sobolev spaces 



Tjn 

HI 

Tjn 



=ff n+3 (T 2 °) x H n+2 (T® ® LG) x tf n+2 (R M ), 
=# n+1 (T 2 °) x H n+1 (T? ® LG) x tf n+1 (R M ), 
=Hl x HI 
=H n {T^) x tf n (LG). 



(3.16) 
(3.17) 
(3.18) 
(3.19) 



One may now show that the operator F in the GRTSE, = 0, defines a 

C 1 -map 

(3.20) 



F : HQ x H™ +2 



H£ for n > 2 



{ 2 } This seems to have been overlooked in [BF]. 
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on the domain of fields which satisfy Condition l.^ 3 ^ For this we need to 

impose suitable but very mild regularity conditions on the unspecified function W 
in (2.5). The linear map L is the first derivative of F in wrt. the second argument: 
D 2 F[$,£\. Ellipticity implies that Image(L) := L(H^ +2 ) C Hg is closed and hence 
H£ splits as orthogonal sum of closed subspaces, given by L's image and the kernel 
of L's adjoint. Hence, since L is self-adjoint, H£ = image(L) © kernel(L). We now 
get an implicit function theorem for the map F if D 2 F[^>, £] is a linear isomorphism, 
i.e. if kernel(L) = {0}. But since L is elliptic, any non-trivial element in the kernel 
may be represented by a C°° function £' which must then satisfy (3.6). Hence a 
trivial kernel is equivalent to the following condition for smooth functions: 

Condition 2. 

fe = kF implies f = 0, k = 0. (3.21) 

Hence we arrive at the following formulation of an implicit-function-theorem for the 
generalized thin-sandwich problem. It may be seen as the analog or generalization 
of Theorem 2 (and 3) in [BF]. 

Theorem. Let n > 2 and f ) G H$ x iL™ +2 be a solution to the GRTSE, 
, f ] = 0, which satisfies Condition 1 (i.e. (2.20-21 ) ) and Condition 2 (i.e. (3.21)). 
Then there exist open neighbourhoods V C HQ of * and W C HQ x H£ +2 of £) 
and a C^map a : V -> Jf™+ 2 such that F[tf',f] = for (#',£') G If <^> G 7 
and f = (T^- 
Consider the action T of R := R — {0} on HQ x H™ +2 , given by T 5 ($, f ) := 
($, 5$, <5£) for 5 G R. It leaves individually invariant the three subsets of points 
(#,£) which 1.) obey Condition 1, 2.) obey Condition 2, 3.) solve the GRTSE. To 
see this, recall that /t is linear in £, hence T$T = 5T. Invariance of the first set 
is now obvious. Further, if fa = kT has only the trivial solution, then so does 
= KT$r, since otherwise (5 -1 ^', k) would be a non-trivial solution to the first 
equation. Hence the second set is invariant. Finally, since V scales with 5 and the 
square-root of expression (2.19) for a with \8\, the GRTSE (2.16) changes at most 
by an overall sign, which proves invariance of the third set. 

The Sobolev embedding theorem for 3-dimensional domains and L 2 -norms implies a continuous 
embedding H n {V) <— > C k (V) for k < n — 3/2. n > 2 is needed to guarantee continuity of the 
functions and gain pointwise control, which is needed in the proof for F being C 1 . 
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We can now repeat the Theorem for each point Tg(^S,^) on the R-orbit of 
(^,0 with open sets Vg, W$ and solution maps a§. In this way the solution map 
a extends to a solution map a * : V* -> Jf™+ 2 , where V* := UaeR^' which 
uniquely represents all solutions in W* := U^gr,^- By construction it satisfies 
cr*($, 8$>) = 5cr($, $) V5 G R. Dropping the superscript *, we formulate this as 

Corollary 1. Let (\I/, £) be as in the Theorem. Then there exist open neighbour- 
hoods V C H$ of U 5e R T <5* and WcF|x Jf ? n+2 of LUr^M*. and a C^-map 
a : V -> i/ c n+2 such that = for (#',£') e W G V and f = <r(V'). 

Moreover, 

a{®,5®) = 5a(®,®), V5 G R (3.22) 

Finally we prove that Condition 2 not only ensures local but also global unique- 
ness. The analogous result has been proven for the traditional RTSE in [BO]/ 4 } 

Proposition 2. Let f ) and £ ) satisfy Condition 1 and the GRTSE. 
Then satisfies the GRTSE <^ there exists a positive function r : E — > R + 

such that r = $ + and f 1 = 6 + /= are related by 

f = rT. (3.23) 

Proof. <=: This follows trivially from the fact that the GRTSE, i.e. equations 
(2.17-18), contain £ only through the combination ^T. 

=^: For s G [0,1], consider the convex combinations £(s) := s£ + (1 — s)£ and 
T s := $ + f^( s ) = sr + (1 — In the following it is useful to think of each T s as 
section in the pulled-back bundle $*T(0) whose fibre at p G £ is a Minkowski space 
pi,5+3iv+M wit}l me t r j c G&fy). Condition 1 requires T(p) and f(p) to be "timelike", 
so that s h- > r s (p) is the straight path connecting these two "timelike" vectors. First 
we show that T(p) and T(p) lie in the interior of the same "light-cone" for some, and 
hence all, p G S. To see this, we consider, for each p, the inner product G(T — f , V) 
with the timelike vector V := g a bg^ of constant length-squared G(V, V) = —3. 
Now, 

/ G(r - f , V) dn = 2 / V a (/T - ^ a ) d/x = 0, (3.24) 

{ 4 } The full statement and proof given in [BO] contains an additional part which is erroneous, as 
was first pointed out in [BF] . If transcribed to our setting, the incorrect part would amount to the 
claim that (3.23) implied r = 1. 
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so that, because E is connected, there exists a point pgS where G(F — f , V)(p) = 0. 
Hence F(p) and F(p) point in the same half of the "light-cone", and so does F s (p), 
since the interior of the half "light-cone" is a convex set. By continuity this must 
then be true at each point peSso that G(T a , F s ) is a negative- valued function on 
S for each s. 

Next we consider the function 

I(s):=S[y,Z(s)}= [ y/-2UG(F s ,F s )d/j. (3.25) 
Jt, 

We have I'(0) = = 7'(1), where ' = £, since £ and f solve the GRTSE. Further- 
more, a straightforward calculation yields: 

I"(s) = [ [2t/] * G(F, F)G(F ± , fx) cfc < 0, (3.26) 



with 



The inequality in (3.26) results from F being "timelike" and F± being "spacelike" 
or zero. But J'(O) = 7'(1) = and 7"(s) < imply I" = 0. On the other hand, 
equality in (3.26) can only be achieved for F± = which is equivalent to (3.23), 
where r must be positive-valued since F and F point in the same half of the "light- 
cone" . □ 

Now (3.23) implies (3.6) with f = f - f and k = so that Condition 2 will 
enforce r = 1 and £ = £. Hence we have 

Corollary 2. 7f f) and £) satisfy the GRTSE and Conditions 1 and 2, then 
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Discussion 



It is obvious that the strategy of the (generalized) thin-sandwich approach cannot 
work for all data. Obvious bad data are those for which the Hamiltonian constraint 
cannot be solved for a nowhere vanishing a. For example, consider fields $ such 
that U > and velocities $ whose gravitational part is pure gauge: g a b = 2V( a £' b ). 
The Hamiltonian constraint implies (V a ?y a ) 2 > Ua 2 , where rj = £' — £, showing 
that a must vanish somewhere since J E d\i V a V a = and E is connected. To avoid 
such situations, Condition 1 or its reversed version, U < and G(T, T) > may be 
imposed. However, if the second condition is chosen formula (3.14) together with 
the proof of Lemma 3 show that L manifestly fails to be elliptic, thus leaving only 
Condition 1. 

The technical Condition 2 has an interpretation in terms of the "canonical data" 
($, V), where aV = $ + /c and where we assume ($, £) to satisfy Condition 1 in 
order to have a^O. Namely, if fa = kT for some non-zero (£', k), then = anV 
says that the same canonical data admit a representation in terms of the new lapse 
function a new = Ka (now possibly with zeros), gauge functions £ new = £' and 
coordinate-velocities $ new = 0. Conversely, if an a new exists such that ct n ew^ = 
/^ now , then f^r = kT with k = a new /a and = ^ n ew Hence Condition 2 precisely 
excludes the existence of other representations of the same canonical data with 
vanishing coordinate-velocities $. We note that Condition 2 may itself be implied 
by simple geometric conditions on $. One such set of conditions is provided by the 
following 

Proposition 3. Condition 2 is implied by Condition 1 and the following conditions 
on $: 



(i) Ric < (Ric = Ricci-tensor of g) (4.1) 



(ii) VaA^ = and A^0 Q = imply = (4.2) 



Proof. Let fe = kT, then G(T, T) < G(fy,fe) < =>• 

J d/x [G abcd 4V {a P' b) V ic P' d) ] =f dfi2 [V [a [3' b] V^(5' h] - R ah [5' a [5 ,b ] < (4.3) 

P' = 0. Since G is positive definite on /e/'s for which (3' = 0, G(f^,f^) < 
implies = which for f3' = is equivalent to the first two equations in (4.2). □ 
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We recall that metrics with Ric < exist on any 3-manifold £ [GY] (e.g. quite in 
contrast to Ric > 0, which is well known to imply a finite fundamental group). (4.2) 
should be read as a mild genericity-condition for the matter fields. For example, 
if we have a single U(l) gauge field and a charged scalar field (here represented 
by a real doublet ((f) 1 7 (p 2 )) then condition (4.2) is satisfied if the scalar field is not 
identically zero. 

Finally we comment on the functional iS[\l/, £] defined in (2.23). Given that 
Condition 1 is satisfied, solutions to the GRTSE are stationary points with respect 
to variations in £ (Lemma 1). Lemma 2 asserts that these must be minima which 
are stable if Condition 2 is satisfied. We now assume Conditions 1 and 2 to hold 
and evaluate S 1 ^, £] at a solution £ = cr(^). We get a C 1 -function S* : HQ — > R + , 

S* [d>, $] = ^ d/i ^^UG^ + f^^ + f^) , (4.4) 

which satisfies 

S*[5V] = \6\S*[V] (4.5) 
for all 5 e R. Standard consequences are 

D 2 S* 8<$>] ($) = sign(5) 5* [$,$], (4.6) 
[$,£$]($,$) = 0, (4.7) 

where in (4.7) we assumed C 2 -smoothness of a. It is tempting to try and regard S* 
as a kind of metric on at least an open subset of the tangent bundle of the space 
of fields. For pure gravity it generalizes a previously considered expression which 
is valid only for constant lapse function [G] and also gives rigorous meaning to a 
formal definition of a distance function given in [CF] . But presently it is unclear to 
us whether S* indeed defines an interesting geometric structure. ^ 

Acknowledgements: I thank Julian Barbour, Robert Bartnik, Gyula Fodor and 
Karel Kuchaf for suggestions and clarifying remarks. This work was financially 

i 5 > One may wonder whether it defined a Finsler metric. For this one would have to show that 
the bilinear form D^S* [<t>, $] is (weakly) non-degenerate. But this is not even the case in finite 
dimensions for functions of the form (4.4) (i.e. sum of square-roots rather than square-root of sum). 

Take e.g. the function f(yi, . . . , y n ) = ■ y/yf- Then didjf = 2difdjf, which is obviously just 

of rank one. 
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